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ABSTRACT  r 

''"■V 

Northern  Hemisphere  kinetic  energy  during  the  FGGE  year  is 

calculated  on  a  twice  daily  basis  at  0000  GMT  and  1200  GMT. 

k 

Graphs  of  time  series  of  total  (K) ,  zonal  (IcJ)  and  eddy  (Kg)  kine¬ 
tic  energy  exhibit  a  strong  annual  cycle  with  superimposed  short 
term  fluctuations  of  about  14  to  16  days  in  K  and  Kg.  This  is  con¬ 
firmed  by  spectral  analysis.  Division  of  Kg  into  standing  and 
transient  components  by  month  shows  good  agreement  with  values  from 
Peixoto  and  Oort  (1974)  for  standing  Kg  but  poorer  agreement  for 
transient  Kg.  A  large  portion  of  the  difference  is  probably  due 
to  interpolation  methods.  Standing  and  transient  Kg  by  month 
and  latitude  show  standing  maxima  associated  with  persistent 

trough  and  ridge  features,  mean  jet  stream  and  summer  monsoon.  A 

.  ri  ■ .  ) 

single  transient  maximum  located  about  5  to  10*  north  of  the  jet- 
relsted  standing  maxlmus  appears  associated  with  cyclone  activity. 
Geographical  plots  of  grid  point  K  showed  Aslan  and  North  American 
east  coast  winter  maxima  near  equal  intensity.  Lower  than  average 
North  American  temperatures  enhanced  barocllnielty,  contributing 
to  the  stronger  than  average  maximum.  Spring  reveals  a  collapsed 
North  American  maxima,  and  the  signature  of  the  Indian  monsoon  ap¬ 
pears  prominent  in  summer.  Available  potential  energy  (A)  quanti¬ 
ties  from  Min  and  Horn  (1982)  allow  K  to  A  ratios  to  be  compared 
with  those  of  Peixoto  and  Oort  (1974). 
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1.  INTRODUCTION 

The  sun  is  the  primary  source  of  energy  to  the  earth  and  its 
atmosphere.  Within  the  atmosphere  energy  appears  in  various  forms 
with  the  transformations  from  one  form  to  another  of  fundamental  sign- 
nificance  to  atmospheric  processes.  Of  the  various  forms,  kinetic 
energy  is  of  most  direct  importance  to  the  meteorologist.  As  the 
energy  of  motion,  it  is  most  directly  related  to  weather  producing 
processes.  Small  vertical  motions  lead  to  cloud  formation  or  dissi¬ 
pation.  Much  larger  horizontal  motions  provide  the  transport  neces¬ 
sary  to  maintain  the  global  energy  balance.  As  a  scalar  quantity, 
kinetic  energy  cannot  describe  motion.  Nevertheless,  a  knowledge 
of  it  is  essential  in  understanding  the  atmosphere's  behavior,  and 
in  particular,  the  atmospheric  energy  cycle.  In  this  study  kinetic 
energy  during  the  year  of  the  First  CARP  Global  Experiment  (FGGE) 
will  be  examined. 

Over  the  years  a  great  deal  of  research  has  been  devoted  to 
kinetic  energy,  its  budget,  and  its  place  in  the  atmospheric  energy 
cycle.  For  the  most  part,  the  research  has  fallen  into  three  broad 
categories:  1)  kinetic  energy  budget  studies  over  limited  areas, 

2)  studies  in  the  wave  number  domain,  and  3)  kinetic  energy  on  a 
hemispheric  scale  as  part  of  the  atmospheric  energy  cycle.  Smith 
(1980)  provides  an  excellent  overview  of  papers  in  the  first  area. 
Saltzmann  (1970)  reviews  earlier  work  in  the  wave  number  domain  and 
Tomatsu  (1979)  lists  more  recent  research.  Oort  (1964)  provides  a 


2 


review  of  kinetic  energy  on  a  global/hemispheric  scale  in  connection 
with  the  atmospheric  energy  cycle.  Peixoto  and  Oort  (1974)  and  Oort 
and  Peixoto  (1974,  1976)  have  performed  extensive  calculations  of 
energy  generation,  reservoirs,  conversions  and  dissipation  using 
data  covering  a  five  year  period.  Although  the  work  here  does  not 
study  the  atmospheric  energy  cycle,  it  is  most  closely  allied  with 
research  in  this  third  category.  Thus,  the  Oort  and  Peixoto  papers 
provide  material  for  comparison  with  kinetic  energy  values  calculated 
here. 

This  thesis  will  examine  Northern  Hemisphere  kinetic  energy 
(K)  during  the  FGGE  year  (December  1,  1978-November  30,  1979)  in  + 

considerable  detail.  Calculations  are  confined  to  the  K  of  hori¬ 
zontal  motions  since  they  are  at  least  four  orders  of  magnitude 
larger  than  the  K  of  vertical  motions.*  The  study  employs  data  in 
isentropic  coordinates  so  as  to  be  compatible  with  the  work  of  Min 
(1981)  and  Min  and  Horn  (1982),  who  compared  available  potential 
energy  (A)  values  computed  using  the  exact  formulation  in  isentropic 
coordinates  with  the  so-called  approximate  form  developed  in  pres¬ 
sure  coordinates.  The  Northern  Hemisphere  K  is  examined  with  respect 
to  both  its  seasonal  and  geographical  distributions.  The  total 
kinetic  energy  (K)  is  partitioned  into  zonal  kinetic  energy  (K^)  and 
eddy  kinetic  energy  (Kg),  with  the  Ke  divided  into  standing  and 

*In  severe  thunderstorms  the  vertical  contribution  might  be  nearly 
as  large,  but  the  larger  scale  focus  and  grid  resolution  employed 
here  minimizes  these  contributions. 
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transient  components.  Time  series  analysis  is  employed  to  gain  in¬ 
sight  into  the  periodicity  of  the  various  K  components.  Using  cal¬ 
culations  from  Min  and  Horn  (1982) ,  ratios  of  K  to  A  are  compared  for 
the  exact  and  approximate  forms  of  A.  Cross-spectral  analysis  of  K 
and  A  are  also  done  to  gain  some  Insight  into  the  relationships  be¬ 
tween  these  energy  forms  during  the  FGGE  year. 

Chapters  2  and  3  contain  a  discussion  of  the  basic  data  and 
calculation  procedures,  respectively.  Following  this.  Chapter  A 
presents  the  seasonal  and  geographical  distributions  of  K  and  ratios 
of  K  to  A.  Chapter  5  presents. the  results  of  time  ser.ies  analysis. 
Chapter  6  contains  a  summary  and  suggestions  for  future  work. 
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2.  BASIC  DATA 

Unlike  many  previous  studies  of  kinetic  energy,  the  FGGE  year 
data  does  not  rely  on  observations  from  only  a  limited  number  of 
radiosonde  stations.  Several  observing  systems  combine  to  provide 
the  best  data  set  currently  available  for  the  entire  globe.  The 
FGGE  data  has  been  and  will  continue  to  be  used  by  many  researchers. 
The  use  of  this  common  data  more  readily  allows  for  comparisons 
between  various  studies.  In  particular,  since  comparisons  are 
planned  between  kinetic  energy  calculations  of  this  study  and  the 
available  potential  energy  results  from  Min  and  Horn  (1982)  and  Min 
(1981),  the  choice  of  the  same  data  period  is  crucial. 

First  consider  a  brief  history  of  the  basic  data  used  here. 

The  National  Meteorological  Center  (NMC)  has  been  generating  global 
data  sets  on  a  twice  daily  basis  for  a  number  of  years.  Townsend 
(1980)  found  that  prior  to  September,  1978  the  gridded  analyses 
generated  by  NMC's  Hough  spectral  analysis  procedure  (Flattery,  1970) 
contained  numerous  cases  of  superadlabatlc  layers.  The  wind  analyses, 
as  expected,  gave  r ona 1 ly- aver aged  meridional  wind  components  that 
were  nearly  xero.  This  effectively  eliminated  the  Hadley  and  Ferrel 
meridional  circulations  he  hoped  to  find.  Thus,  Townsend  was  unable 
to  use  the  Hough  gridded  data.  In  September,  1978,  in  preparation 
for  j  FGGE  year,  NMC  introduced  an  optima  interpolation  global 
data  assimilation  system  (Bergman,  1979;  McPherson  et  al.,  1979). 
Townsend’s  subsequent  testing  showed  only  an  occasional  superadia- 
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bade  layer*  and  Che  wind  field  now  contained  non-zero  meridional 
wind  components.  Grldded  data  based  on  this  scheme  were  used  by 
Townsend  (1980)  and  Min  (1981)  and  will  be  employed  In  this  study. 

a.  The  FGGE  level  Ilia  data  set 

The  FGGE  level  Ilia  data  set  used  In  this  study  was  obtained 
from  the  National  Center  for  Atmospheric  Research  (NCAR).  The  data 
Include  pressure,  temperature,  wind  and  relative  humidity  compiled 
twice  dally  at  0000  GMT  and  1200  GMT.  The  data  were  arrayed  In  a 
2.5*  latitude-longitude  grid  for  each  of  the  12  mandatory  isobaric 
levels:  1000,  850,  700,  500,  400,  300,  250,  200,  150,  100,  70,  and 
50  mb.  A  heterogeneous  collection  of  observing  systems  provided  the 
raw  Information.  The  systems  included  conventional  radiosonde  ob¬ 
servations  (BAOBs),  surface  reports  from  land  stations  and  ships 
and  conventional  aircraft  reports  (AIREPs).  A  unique  feature  of 
FGGE  was  the  extensive  coverage  of  data  sparse  oceanic  regions  by 
satellite.  T>e  geostationary  satellite  system  generated  cloud- 
tracked  winds  and  a  polar  orbiting  satellite,  TIR0S-N,  yielded  verti¬ 
cal  temperature  profiles.  In  tropical  regions,  two  special  systems 
provided  observations:  the  Tropical  Wind  Observing  Ships  (TWOS)  and 
Tropical  Constant  Level  Balloons  (TCLB).  During  the  two  FGGE  Spe¬ 
cial  Observing  Periods  of  January-February  and  May- June,  OMEGA- 
systam  dropsondes  produced  temperature,  relative  humidity  and  wind 
observations  along  specific  aircraft  routes  primarily  in  the  tropics. 
See  Fleming  et  al.  (1979a, b)  for  a  more  complete  description  of  the 
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data  sec. 

The  variety  of  observing  systems  produced  meteorological 
measurements  with  differing  error  characteristics.  Consequently, 

NMC  developed  a  data  assimilation  system  based  on  a  statistical 
interpolation  scheme.  The  statistical  method,  an  optimum  interpo¬ 
lation  (01)  analysis  scheme,  is  capable  of  multivariate  analysis 
of  winds  and  temperatures.  The  01  scheme,  which  assumes  that  both 
observed  and  forecast  values  are  subject  to  error,  updates  pre¬ 
dicted  values  with  observations.  For  the  details  of  the  NMC  global 
assimilation  system  see  Bergman  (1979)  and  McPherson  et  al.  (1979). 

b.  Isobaric-isentropic  data  processing 

In  processing  the  data,  NMC  generated  the  global  analysis 
in  sigma  coordinates  and  then  used  vertical  interpolation  from  sigma 
to  lsobaric  coordinates  to  obtain  the  isobarlc  analysis.  The  pro¬ 
cessed  Northern  Hemisphere  data  were  stored  in  145x37  longitude- 
latitude  arrays.  This  corresponds  to  a  2.5s  longitude-latitude 
grid.  Data  arrays  were  generated  for  the  12  mandatory  pressure 
levels  from  1000  to  50  mb  on  a  twice  daily  basis  at  0000  GMT  and 
1200  GMT. 

Of  a  possible  730  observation  times  during  the  FGGE  year,  703 
were  available  in  the  data  set.  Min  (1981)  eliminated  two  addi¬ 
tional  periods  because  of  suspect  data  in  his  available  potential 
energy  study.  While  calculating  kinetic  energy,  two  more  time 
periods  were  eliminated  because  of  unrealistic  wind  values  (in  excess 
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of  1000ms  *")  at  several  grid  points.  Thus  a  total  of  699  observa¬ 
tion  times  were  employed.  Because  the  missing  periods  shown  In  Table 
1  were  about  equally  distributed  throughout  the  12  months,  no  attempt 
was  made  to  Interpolate  for  missing  analyses. 

Min  (1981)  made  some  Initial  calculations  of  available  poten¬ 
tial  energy  to  test  different  horizontal  and  vertical  resolutions. 
From  his  calculations  It  was  decided  to  employ  a  5°  latitude- 
longitude  grid  and  10K  lsentroplc  vertical  intervals.  This  study 
employs  the  same  horizontal  and  vertical  resolution. 

Obtaining  lsentroplc  data  from  the  NMC  analyses  proved  to  be  a 

major  undertaking .  Townsend  (1980)  generated  the  lsentroplc  data 

by  vertical  interpolation  of  mandatory  level  data  to  lsentroplc 

surfaces.  He  used  an  algorithm  developed  by  T.  Whittaker  at  the 

University  of  Wisconsin.  The  algorithm  assumes  a  linear  variation 

1/k 

of  pressure  while  all  other  parameters  vary  with  8  .  The  ten 

mandatory  pressure  levels  from  1000  to  100  mb  were  used  in  this 
interpolation.  The  top  lsentroplc  surface  is  at  370K,  which  corres¬ 
ponds  roughly  to  the  150 mb  level,  and  thus,  includes  essentially  all 
of  the  tropospheric  circulation.  Potential  temperature  (6)  extends 
downward  from  370K  in  10K  increments.  The  lowest  6  surface  present 
during  a  data  period  depends  on  the  minimum  observed  potential 
temperature  at  the  earth's  surface  at  any  grid  point  over  the  hemis¬ 
phere.  Values  as  low  as  220K  were  observed  during  the  winter.  Thus, 
data  from  the  ten  mandatory  pressure  levels  were  interpolated  to  as 
many  as  16  lsentroplc  surfaces  extending  from  370K  to  as  low  as  220K. 
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Table  1.  Distribution  of 

Ilia  data  set. 

the  number 

of  missing  data 

in  FGGE  level 

0000 

GMT 

1200 

GMT 

Total 

December,  1978 

1 

0  (2)* 

1  (3) 

January,  1979 

1 

1 

2 

February 

1 

0 

1 

March 

2 

1 

3 

April 

1 

1 

2 

May 

0 

2 

2 

June 

1 

3 

4 

July 

0 

1 

1 

August 

1 

1 

2 

September 

4  (6)* 

3 

7  (9) 

October 

1 

0 

1 

November,  1979 

1 

0 

1 

Total 

14  (16) 

13  (15) 

27  (31) 

*Two  time  periods  in  December  were  eliminated  because  of 
excessive  winds  at  some  grid  points.  The  two  additional  periods 
eliminated  in  September  resulted  from  very  small  A  values  found  by 
Min  (1981). 


3.  EQUATIONS  AND  COMPUTATIONAL  METHODS 
a.  Kinetic  energy  equations 

In  this  study,  total  kinetic  energy  (K)  per  unit  area  Is  de¬ 
fined  as 

11  •  ife  Js  J~9]  <“2  +  y2)  T*  i,liS  (1) 

where  g  is  gravitational  acceleration,  S  is  hemispheric  surface 
area,  u  and  v  are,  respectively,  the  eastward  and  northward  compon¬ 
ents  of  the  wind,  p  is  pressure,  6  is  potential  temperature  and 
2 

dS  -  a  coS(j)d<j>dX,  where  a  is  the  earth's  radius,  $  is  latitude  and 
X  is  longitude.  The  vertical  integration  performed  in  isentroplc 
coordinates  extends  from  220K  to  370K. 

K  was  partitioned  into  zonal  (Kg)  and  eddy  (Kg)  components: 

*z  ■  ife  JSL;  «“i2 + w2>  n dMs  <2> 


h 


(u*2  +  V*2)  d6dS 


(3) 


where  [  ]  indicates  a  zonal  average  and  ()  is  the  deviation  from  the 
zonal  average.  Appendix  A  contains  detailed  information  on  the 
actual  computational  procedures. 
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Ic  is  instructive  to  further  divide  into  its  standing  and 
transient  components.  The  transient  component  of  Kg  (ICEtr)  may  be 
written  as 

*Etr  “  *E  “  *Est 

where  the  overbar  implies  a  time  average.  The  subscripts  tr  and  st 
identify  the  transient  and  standing  components.  is  given  by 

*E.t  •  i  JJ  *  <«*  ♦  '*  >  H dMS  (5) 

S  0s 

In  the  above  equation,  isobaric  data  is  first  time  averaged  for  each 
ronth  and  then  interpolated  to  0  coordinates.  This  method  is  em¬ 
ployed,  rather  than  time  averaging  the  twice  daily  isentropic  data, 

» 

since  the  height  of  a  pressure  surface  fluctuates  less  than  the 
height  of  an  isentropic  surface  during  a  month. 

The  geographical  distribution  of  K  and  its  standard  deviation 
is  also  investigated.  Each  5*  latitude-longitude  box  can  be  repre¬ 
sented  by  its  vertically  integrated  grid  point  value.  Specific  units 
(jtg-^)  are  employed  because,  when  calculating  the  vertical  Integral 
in  9  coordinates,  the  mass  between  0g  and  0t  can  vary  from  one  point 
to  another.  The  K  equation  in  specific  units  is 

XJ,'  <-2  ♦’*>$**• 

s  qa _ 

r  (p  -  -  p„  )d* 

J  8  ' *s£c  rtop 


K 


(6) 


where  Che  denoalnacor  represents  Che  mass  in  the  column  determined 
from  the  surface  pressure  and  Che  pressure  at  Che  cop  of  che  domain 
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from  the  surface  pressure  and  Che  pressure  at  Che  cop  of  che  domain, 

i.e.,  p  at  370K,  and  small  s  Is  Che  area  of  the  grid  box.  See 

cop  .  , 

Appendix  A  for  a  discussion  of  area  calculation.  The  equations 

for  Kg  and  Kg  in  specific  units  are  similar  to  Equation  6.  Mote 

that  since  che  wind  is  broken  inCo  a  zonal  average  and  a  deviation 

* 

at  each  grid  point  (e.g. ,  u  *  [u]  +  u  ),  a  cross  product  term  is 

2 

implicitly  included  in  grid  point  values  of  K.  For  example,  u  • 

2  *  *2  2  *2 
[u]  +  2[u]  u  +  u  .  Since  Kg  a  [u]  and  Kg  cc  u  ,  the  cross  product 

*  * 
term  2[u]u  makes  neither  a  zonal  nor  eddy  contribution.  The  2[u]u 

term  integrates  to  zero  in  zonal  rings  and  in  the  hemispheric  total; 

however,  it  is  contained  implicitly  in  the  total  K  value  at  each 

grid  point.  Thus,  at  grid  points,  K  is  not  equal  to  the  sum  of  Kg 

and  Kg. 


b.  Time  series  analysis 


A  wealth  of  information  about  a  time  series  can  be  readily  ob¬ 
tained  by  performing  a  Fourier  analysis.  The  Fourier  analysis  fits 
sine  and  cosine  waves  to  the  data  to  estimate  the  contribution  of 
each  Fourier  frequency  to  Che  actual  time  series.  The  variance  ex¬ 
plained  by  a  particular  Fourier  frequency  can  be  obtained  from  the 
power  (or  variance)  spectrum  which  is  computed  as 

1  2  2 

Pi  "  I  (ai  +  bi  >  (7) 


(7) 
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where  each  gives  the  amount  of  variance  explained  by  a  particular 
Fourier  frequency.  By  standardizing  the  data  prior  to  analysis, 
represents  a  percentage  of  total  variance  explained  at  the  i'th 
frequency.  The  coefficients,  and  are  given  by 

2  N 

a  -  ~  Z  (X(t)sin(iu>t)  ] 

1  N  t-1 


I  [X(t)cos(ia)t) ] 

1  N  t-1 

where  1  may  be  any  Integer  from  1  to  y  -  1.  The  summation  is  over 
the  N  observations,  and  X(t)  la  the  value  of  the  aeries  at  point  t. 

The  frequency,  oj,  is  equal  to  2ir/T,  where  T  is  the  period  of  record. 
By  subjecting  two  simultaneous  time  series  to  spectral  analysis, 
one  learns  how  the  two  series  vary  together.  This  cross  spectral 
analysis  consists  of  two  parts,  the  coincident  spectrum  (C)  and  the 
quadrature  spectrum  (Q).  C  measures  the  amount  of  the  two  series  in- 
phase  relation  and  Q  measures  the  out-of-phase  relation.  The  coin¬ 
cident  and  quadrature  spectrums  are  given  by 


C 

Q 


i 

i 


a. a  .  +  b  .  b  , 
xi  yi  xi  yi 


b  ,  a  .  -  a  -  b  _ 
xi  yi  xi  yi 


(9) 


where  the  subscripts  denote  the  coefficients  from  the  two  simultan¬ 
eous  series  X  and  Y  for  the  i'th  frequency.  A  phase  angle  differ¬ 
ence  (A8)  is  computed  from  C  and  Q  at  each  frequency  as 


tan(AfJj) 


(10) 


Finally,  to  determine  how  good  a  relationship  exists  between 
the  two  series,  the  coherence  (CH)  is  computed  as 


The  coherence  is  estimated  by  summing  over  a  finite  band  of  fre¬ 
quencies.  This  is  frequently  done  for  estimating  P  as  well  since 
this  increases  the  degrees  of  freedom  for  the  spectral  estimate. 

This  will  be  discussed  more  completely  in  Chapter  5.  Coherence  can 
range  from  0  to  1  and  is  analogous  to  the  square  of  a  correlation 
coefficient  except  that  CH  is  a  function  of  frequency. 

A  more  complete  discussion  of  spectral  analysis  can  be  found 
in  Panofsky  and  Brier  (1958)  and  Bendat  and  Plersol  (1971).  Addi¬ 
tionally,  Appendix  B  contains  a  discussion  of  the  procedures  employed 
prior  to  using  a  fast  Fourier  transform  (FFT)  to  obtain  the  spectral 
results  presented  in  Chapter  5. 


t 
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4.  RESULTS 

This  section  presents  the  findings  derived  from  the  tvice  daily 
kinetic  energy  calculations  performed  for  the  Northern  Hemisphere 
using  the  level  Ilia  FGGE  data  set.  The  time  series  of  kinetic 
energy  values  are  discussed  first,  followed  by  a  more  detailed  look 
at  Kg  which  is  partitioned  into  its  standing  and  transient  compo¬ 
nents.  Next,  the  spatial  distribution  of  the  mean  seasonal  K  and 
its  standard  deviations  are  reviewed.  Finally,  ratios  are  calcu¬ 
lated  using  A  from  Min  and  Horn  (1982) . 

a.  Time  series  of  kinetic  energy  values 

Figure  1  shows  the  twice  daily  values  of  the  total  kinetic 

energy  (K)  and  its  zonal  (Kg)  and  eddy  (Kg)  components  vertically 

integrated  from  220K  to  370K  over  the  area  of  the  Northern  Hemi- 

5  -2 

sphere.  Units  are  xlO  Jm  .  All  three  curves  exhibit  a  strong 
annual  cycle  which  responds  to  the  solar  cycle,  with  the  zonal  term 
showing  slightly  more  seasonal  variation  than  the  eddy  term.  This 
is  in  agreement  with  the  findings  of  Peixoto  and  Oort  (1974)  and 
is  what  one  expects  from  examining  upper  air  charts  over  the  course 
of  a  year.  K  is  highest  during  the  winter  months  when  radiation 
differences  are  largest  between  the  equator  and  pole.  This  leads  to 
large  north-south  temperature  contrasts,  which  yield  strong  zonal 
winds  aloft  via  the  thermal  wind  relationship.  Excessively  large 
thermal  gradients  and  zonal  winds  would  result  were  it  not  for  the 
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Figure  1.  Time  series  of  twice  daily  vertically  Integrated  K  for 
the  Northern  Hemisphere  during  the  FGGE  year.  Total 
(dash-dot) ,  zonal  (dash)  and  eddy  (solid).  Units: 

105  Jm-2. 


development  of  vigorous  eddy  motions  which  transport  heat  poleward. 

The  K  curve  exhibits  a  strong  pulsing  in  the  winter  on  the  order  of 
12-16  days.  This  strong  pulsing  returns  in  late  autumn  when  radia¬ 
tion  deficits  in  high  latitudes  are  Increasing.  Visual  inspection 
shows  that  when  K  exhibits  peaks,  so  does  K^;  thus  eddy  motions  appear 
to  be  the  primary  contributor  to  the  short  term  shape  of  the  K 
curve,  while  zonal  motions  tend  to  provide  a  quasi-steady  background 
state. 

As  spring  approaches,  a  nearly  steady  decline  in  K  begins  and 
extends  through  mid-summer .  The  decrease  is  most  evident  in  V 
This  is  a  reflection  of  the  decreasing  baroclinicity  between  the 
tropics  and  high  latitudes.  The  decline  in  Kg  is  less  pronounced. 
East-west  heating  differences,  which  result  in  eddy  motions,  are 
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evidently  becoming  relatively  more  Important  and  partly  replace 
the  vigorous  winter-time  eddies  which  are  required  for  a  poleward 
energy  transport.  Consequently,  Kg  decreases  less  than  Kg  during 
the  warmer  seasons. 

Summer  is  a  period  of  low  K  values,  but  as  autumn  approaches, 

K  begins  to  Increase.  The  short  period  fluctuations  return  in  K 
and  Kg  while  Kg  slowly  Increases.  Mote  that  the  FGGE  year  ended 
at  a  somewhat  lower  kinetic  energy  level  than  it  began.  The  dif- 
ference  is  made  up  almost  entirely  by  lower  levels  of  Kg  at  year's 
end.  •  • 

Monthly  mean  values  of  twice  dally  kinetic  energy  are  shown 
in  Figure  2.  Monthly  values  are  connected  to  show  the  overall  trend 
from  month  to  month.  Note  that  the  FGGE  year  begins  with  Kg  and  KE 
about  equal  and  ends  (November)  with  Kg  significantly  greater  than 
Kg.  Figure  3  is  identical  to  Figure  1  except  that  it  shows  the 
time  series  of  0000  GMT  and  1200  GMT  kinetic  energy  separately. 
Although  not  evident  from  the  figures,  diurnal  changes  in  K  are 
small.  Unlike  the  hemispheric  curves  of  Min  and  Horn  (1982)  and 
the  regional  findings  of  Koehler  and  Min  (1983)  for  A,  which  is 
strongly  influenced  by  the  variability  of  temperature  near  the 
surface,  there  is  no  large  dally  variation  in  wind  that  contributes 
strongly  to  the  vertically  integrated  K. 


Figure  2. 


Monthly  means  of  K  values  shown  in  Fig.  1.  Total  (dash 
dot),  zonal  (dash)  and  eddy  (solid).  Units:  105  JnT2. 


Figure  3.  Tima  series  of  dally  K  values  from  Fig.  1  but  for  only 
a)  0000  GMT  end  b)  1200  GMT.  Total  (dash-dot),  zonal 
(dash)  and  eddy  (solid).  Units  10^  Jm~2. 
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b.  Standing  and  transient  components  of  Kg 

Special  attention  is  given  to  Kg  since  it  is  the  energy  of 
motion  in  atmospheric  disturbances.  Kg  was  partitioned  into  its 
standing  (Kgst)  and  transient  (Rgcr)  components  using  Equations 
4  and  5.  Figure  4  shows  the  monthly  values  of  Kgat  and  Kg  . 

shows  only  a  weak  annual  cycle  and  compares  favorably  with 
the  findings  of  Peixoto  and  Oort  (1974) .  The  Kgtr  curve  is  larger 
than  Kgst  throughout  the  year  and  exhibits  a  stronger  annual  cycle. 
The  Kgtr  values  found  here  are  about  27Z  lower  than  those  of  Peixoto 
and  Oort.  From  examination  of  Table  2  it  is  clear  that  the  majority 
of  the  monthly  differences  in  total  kinetic  energy  (K)  between  the 
studies  seems  to  be  accounted  for  by  the  differences  in  Kgtr>  The 
values  in  Table  2  are  in  specific  units  (Jkg  S ;  for  comparison  pur¬ 
poses  Peixoto  and  Oort  (1974)  values  were  converted  to  specific 
units.  This  was  done  because  of  the  different  vertical  domains, 
l.e.,  1000-100  mb  for  Peixoto  and  Oort  and  surface  to  370K  (about 


d  ■  i - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1  b 


°oec  Fta  ofn  mfr  jun  jul  5uo  s?p  oct  nov 

Figure  4.  Eddy  K  (dash-dot)  partitioned  into  standing  eddy  (dash) 
and  transient  eddy  (solid).  Units  10^  Jm~*. 
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Table  2.  Monthly  mean  values  of  kinetic  energy.  Results  of  this 

study  are  compared  with  the  five  year  average  values  from 
Peixoto  and  Oort  (1974)  (converted  to  specific  units) . 
Units :  Jkg"^- 


K 

Peixoto 

*Z 

Peixoto 

Mo. 

FGGE 

&  Oort 

Mo. 

FGGE 

&  Oort 

D 

147.2 

167.7 

D 

72.7 

73.7 

J 

158.2 

188.8 

J 

75.9 

84.4 

F 

153.7 

185.1 

F 

77.6 

90.4 

M 

143.9 

166.5 

M 

72.0 

79.7 

A 

111.2 

132.7 

A 

48.4 

54.4 

M 

91.7 

102.4 

M 

33.4 

32.3 

J 

66.4 

84.1 

J 

20.2 

22.2 

J 

59.0 

73.3 

J 

18.9 

20.4 

A 

58.8 

80.3 

A 

18.4 

24.9 

S 

68.0 

91.4 

S 

23.9 

29.4 

0 

92.2 

111.1 

0 

35.9 

37.1 

N 

126.0 

138.0 

N 

47.9 

53.8 

rt 

Peixoto 

*Est 

Peixoto 

Mo. 

FGGE 

&  Oort 

Mo. 

FGGE 

&  Oort 

D 

53.9 

75.0 

D 

20.5 

18.9 

J 

61.1 

78.3 

J 

21.2 

23.1 

F 

57.3 

72.6 

F 

18.9 

22.1 

M 

54.4 

70.8 

M 

17.5 

16.1 

A 

50.0 

67.9 

A 

12.7 

10.5 

M 

43.2 

60.9 

M 

15.1 

9.4 

J 

34.4 

51.1 

J 

11.8 

10.9 

J 

29.5 

44.3 

J 

10.5 

11*7 

A 

30.4 

44.9 

A 

10.1 

10.5 

S 

34.3 

52.2 

S 

9.8 

9.8 

0 

42.6 

62.6 

0 

13.7 

11.3 

N 

56.2 

68.4 

N 

21.8 

15.8 
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150  mb)  here.  Although  the  calculations  obtained  in  this  study  do 

not  Include  contributions  from  about  150  to  100  mb,  Kgtr  still  seems 

lower  than  what  might  be  expected  from  Peixoto  and  Oort's  five  year 

average  data  from  May  1958  to  April  1963.  The  remaining  difference 

can  in  part  be  explained  by  comparing  calculation  methods.  Peixoto 

2 

and  Oort  interpolate  monthly  statistics  (such  as  v'  )  from  stations 
to  grid  points  rather  than  interpolating  the  basic  data  (e.g.,  v) 
and  then  calculating  energy  quantities  as  was  done  here.  The  re¬ 
sult  is  that  the  estimates  are  different. 

To  illustrate  the  calculation  differences,  assume  a  square 
array  of  equally  spaced  data  points  surrounding  a  central  grid 
point.  The  data  points  are  numbered  clockwise  starting  in  the 
upper  left  hand  corner;  the  central  grid  point  is  zero.  As  an 
example,  imagine  a  sinusoidal  wave  propagates  through  this  array 
in  five  time  periods,  a.  .  through  e.,  where  both  time  periods  a. 
and  e.  correspond  with  the  data  grid  positioned  in  the  light  winds 
under  the  ridge.  Period  c.  centers  the  grid  directly  on  the  trough 

while  time  periods  b.  and  d.  place  the  grid  at  intermediate  posi- 

2 

tions.  With  this  arrangement,  v'  is  computed  by  the  two  different 

methods.  The  upper  portion  of  Table  3  shows  the  hypothetical  v 

components  at  the  four  data  points  and  their  deviations  from  the 

time  average  (v')  in  parentheses.  Values  for  the  central  grid 

point  Oare  the  result  of  Interpolating  at  each  time  period.  In 

2 

the  lower  portion  of  the  table,  v'  calculated  from  the  interpolated 
2  -2 

v'  values  is  2A  m  s  .  The  values  to  the  extreme  right  are  the  re- 


comparison  of  two  different  interpolation  techniques. 

fnaolal  AVA  -L _ _  M.  a  _  -  .  .  .  • 
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suit  of  cine  averaging  at  each  data  point  and  then  interpolating. 

2  2  -2 

For  v'  this  value  is  101  m  s  .  As  shown  by  this  hypothetical 
example,  the  method  of  interpolation  gives  large  differences  for  the 
v*  quantity.  The  computational  method  chosen  by  Peixoto  and  Oort 
will  produce  larger  (or  equivalent)  K£tr  estimates.  K£st  values 
are  comparable  between  studies  since  the  interpolation  was  anal¬ 
ogous  to  that  shown  in  the  upper  portion  of  Table  3,  where  no  dif¬ 
ferences  result.  In  addition,  Peixoto  and  Oort  find  Kg  as  the  sum 
of  KggC  and  Kgcr>  Here,  is  found  as  a  residual  of  the  average 

of  twice  daily  values  of  Kg  and  monthly  Kggt  (Eq.  4).  Thus,  the 
two  methods  will  give  dissimilar  values  for  Kgtr  and  K_.  These 
differences  should  be  tacitly  noted  when  comparing  Kg  values. 

To  further  understand  the  behavior  of  Kg>t  and  Kg£r,  plots 
showing  each  as  a  function  of  time  and  latitude  were  prepared. 

Kgst  and  Kgtr  ere  presented  in  Figure  5.  The  average  values  for 
each  zonal  ring  (weighted  by  the  area  of  the  ring)  are  plotted  for 
each  month.  The  standing  component  of  Kg  in  Figure  5a  exhibits 
three  relative  maxima  during  the  year.  These  regions  are  marked 
by  heavy  dashed  lines.  The  first  region  is  a  band  of  strong  Kg#t 
between  55-60 *H  present  during  the  winter  months.  This  maximum 
corresponds  closely  with  the  latitude  of  the  east-vest  monthly  mean 
trough  positions.  The  mesa  /00  mb  height  map  from  Taubensee  (1979) 
reproduced  in  Figure  6  for  December  shows  s  well  developed  low  over 
northeast  Siberia  and  a  trough  extending  eastward  through  a  south¬ 
ward  displaced  Icelandic  low  to  just  west  of  Ireland.  The  east-west 


axis  of  these  features  Is  roughly  along  55“N.  These  features  contri¬ 


bute  to  Kg  by  large  deviations  in  both  u  and  v.  The  high  latitude 
maximum  persists  through  January,  February  and  March  and  exhibits 
behavior  consistent  with  December.  The  maximum  between  70  and  80“ N 
which  appears  in  February  overshadows  the  relative  maximum  near  60“N. 
An  examination  of  the  mean  700  mb  map  from  Dickson  (1979) ,  repro¬ 
duced  in  Figure  7,  shows  a  strong  omega  block  over  the  Bering 
Strait  and  a  wave  number  2  pattern  with  intense  upper  level  lows 
over  north  central  Asia  and  the  Canadian  archipelago.  These  features 
contribute  very  strong  departure  values.  Inspection  of  daily  500  and 
300  mb  maps  from  the  European  Centre  for  Medium  Range  Forecasts 
(Bjorhelm  et  al. ,  1981)  also  reveal  patterns  that  would  contribute 


Figure  6.  Mean  700  mb  height  contours  (dam)  for  December  1978  from 
Taubensee  (1979). 
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Figure  7.  Mean  700  mb  height  contours  (dam)  for  February  1979 
from  Dickson  (1979). 

to  a  large  standing  eddy  component.  From  this  comparison  it  seems 
plausible  that  strong  and  persistent  high  latitude  trough  and  ridge 
positions  are  delineated  by  large  values  of  K£sC>  particularly  when  . 
they  occur  in  essentially  the  same  zonal  ring. 

The  second  region  of  large  K£st  can  be  followed  through  nearly 
the  entire  FGGE  year.  Its  seasonal  cycle  very  closely  follows  the 
latitude  bands  of  the  axis  of  the  mean  jet  stream  position.  This 
agrees  favorably  with  the  position  of  maximum  zonally  averaged 
monthly  geostrophic  winds  computed  from  the  pressure  level  height 
information  from  Crutcher  and  Me serve  (1970).  This  maximum  will  be 
further  discussed  below.  It  should  be  noted  that  ao  double  maximum 
appears  in  November  as  had  been  present  in  the  previous  December. 
Inspection  of  the  700  mb  flow  pattern  presented  by  Dickson  (1980) 


26 

does  not  show  as  broad  or  as  strong  700  mb  features  as  were  present 
early  In  the  FGGE  year. 

the  third  region  of  relatively  high  values  of  Kgst  occurs  in 
summer  at  about  15°N.  This  feature  is  probably  associated  with  the 
land-ocean  heating  contrasts  which  generate  the  summer  monsoon. 

Some  evidence  that  this  is  a  monsoon  signal  will  be  given  in  the 
next  section. 

The  time-latitude  plot  of  transient  Kg  shown  in  Figure  5b  ex¬ 
hibits  only  one  relative  maximum  during  the  year.  The  heavy  dashed 
line  drawn  through  the  figure  appears  similar  to  the  second  maximum 
discussed  for  Kgst.  but  is  more  clearly  defined.  It  tends  to  be 
located  about  5  to  10*  north  of  the  KgsC  maximum  for  each  month 
throughout  the  year.  As  a  possible  explanation  of  this  displace¬ 
ment  one  might  imagine  a  transient  eddy  (e.g. ,  a  cyclone)  forming 
near  the  axis  of  the  Kgat  maxima  (near  a  jet  stream  maxima).  As 
the  storm  matures  it  moves  to  the  northeast  and  intensifies.  This 
process  of  an  eddy  moving  northeast  and  Intensifying  would  tend  to 
move  the  maximum  Kgcr  contribution  to  a  more  northerly  latitude  band. 
A  comparison  of  Kgtr  with  the  latitudes  of  FGGE  year  cyclogenesis 
prepared  by  L.  Whittaker,  and  shown  in  Figure  8,  lends  support  to 
this  view.  Comparison  of  Figures  5b  and  8,  in  general,  shows  that 
the  maximum  values  of  Kgcr  appear  5  to  10*  north  of  the  average 
latitude  of  greatest  cyclogenetlc  activity  as  determined  from  the 
12  cyclone  contour.  Noteable  exceptions  occur  In  February  and  May 
when  the  maxima  cyclogenesis  count  in  one  latitude  ring  occurs  north 
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of  Che  &£Cr  maximum.  Comparison  of  FGGE  cyclogenesis  with  chat 
found  by  Olson  (1982)  for  the  period  1958-1977  shows  cyclogenesis 
during  February  and  May  of  the  FGGE  year  to  be  displaced  about  10° 
north  of  the  long  term  mean.  These  two  months  seem  to  contradict 
the  Idea  that  K^tr  should  appear  north  of  the  cyclone  maximum,  how¬ 
ever,  it  should  be  noted  that  cyclone  Intensity  is  not  represented 
In  Figure  8.  It  is  possible  that  a  lesser  number  of  long  lived. 
Incense  cyclones  in  a  lower  latitude  belt  could  contribute  to  larger 
Kgcr  than  frequent,  short  lived,  weak  cyclones  in  a  higher  latitude 
belt.  An  examination  of  the  cyclones  that  formed  in  high  latitudes 
during  February  and  May  revealed  the  storms  were  generally  shorter 
lived  and  less  intense  than  those  further  south  (L.  Whittaker, 

1983) .  In  addition,  no  mention  has  been  made  regarding  how  anti- 

\ 

cyclones  contribute  to  Kgtr. 

Min  and  Horn  (1982)  showed  that  transient  eddy  A  (Agtr)  in 
the  exact  form  for  the  FGGE  year  corresponded  quite  closely  to  the 
zones  of  maximum  cyclogenesis  found  by  Whittaker  and  Horn  (1981) 
for  North  America  and  Whittaker  and  Horn  (1982)  and  Klein  (1957) 
for  the  Northern  Hemisphere.  It  is  Interesting  to  compare  Kgtr  with 
transient  eddy  A  from  their  study.  Figure  9  shows  the  transient 
eddy  A  from  Min  and  Horn  (1982)  for  the  exact  and  variable  approxi¬ 
mate  formulations  of  A.  Comparing  Figure  9a  with  Figure  5b  reveals 
the  similar  shape  of  the  maxima  of  both  plots.  The  maxima 

appears  at  or  slightly  north  (within  one  5*  latitude  ring  on  average) 
of  the  exact  Agtr  maximum.  If  cyclogenesls  is  closely  associated 
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with  the  A^tr  maximum,  then  once  formed,  the  cyclones  Intensify  and, 
generally,  move  to  the  northeast.  Thus,  they  displace  the 
maximum  northward.  In  contrast,  comparison  of  the  approximate  Agtr 
maximum  (Fig.  9b)  with  the  Kgtr  maximum  shows  a  poor  fit.  In  fact, 
Kgtr  averages  5  to  10s  south  of  the  Ag maximum.  This  seems  to 
further  support  the  argument  of  Min  and  Horn  (1982)  that  the  exact 
formulation  of  A  should  be  used  for  studying  cyclone  activity. 

c.  Geographic  distribution  of  K,  Kg  and  their  standard 
deviations 

In  this  section  the  geographical  distributions  of  K  and  Kg  are 
presented  along  with  their  standard  deviations.  The  values  are  ex¬ 
pressed  in  specific  units  (Jkg  1)  because  the  vertical  integration 
does  not  contain  the  same  amount  of  mass  in  every  grid  square.  Dis¬ 
tributions  from  Lahey  et  al.  (1960)  of  300  mb  geostrophlc  K  are 
referenced  in  this  section  to  provide  a  rough  estimate  of  differences 
in  the  FGGE  year  distributions  from  a  multi-year  average.  Because 
the  integrated  values  of  K  versus  geostrophlc  K  at  one  level  are 
not  entirely  comparable,  only  very  general  comparisons  are  possible. 
However,  it  should  be  noted  that  300  mb  is  near  the  level  of  maximum 
K  contribution.  Mid-season  months  from  Lahey  et  al.  (1960)  were 
chosen  for  comparison  with  seasonal  maps  produced  here. 

Figure  10a  shows  K  for  grid  squares  and  cyclone  tracks  for  the 
FGGE  winter  (DJF) .  Two  strong  centers  of  K  appear,  one  off  the  east 
coast  of  Asia  and  an  equally  strong  one  Just  north  of  Cape  Hattecas 
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on  the  U.S.  coast.  A  third,  weaker  maximum  appears  over  the  Arabian 
peninsula.  These  are  generated,  primarily,  by  the  strong  upper 
level  jet  streams  in  these  regions  of  strong  thermal  land/ocean 
contrast.  The  position  of  the  centers  are  in  general  agreement  with 
Lahey  et  al.  (1960),  however,  the  two  FGGE  east  coast  maxima  are 
nearly  the  same  intensity  while  the  long  term  average  shows  the 
Asian  maximum  to  be  much  stronger  than  the  North  American  maximum. 
Extremely  cold  weather  over  the  North  American  continent  during  the 
FGGE  winter  probably  led  to  larger  baroclinicity  which  supported 
stronger  winds  and  gave  the  larger  K  maximum.  Standard  devia¬ 
tions  of  K  in  Figure  10b  look  similar  to  the  map  of  K.  The  figure 
seems  to  reflect  the  pulsing  and  east-west  elongation  of  the  wind 
through  its  average  position. 

In  Figure  10a,  cyclone  tracks  from  the  FGGE  winter  agree  quite 
well  with  observations  of  cyclone  activity  in  relation  to  upper 
level  winds.  Since  K  is  vertically  integrated,  and  the  largest 
contributions  occur  from  the  layers  near  wind  maxima,  one  can  think 
of  the  K  map  roughly  representing  mean  jet  stream  positions.  The 
superimposed  cyclone  tracks,  in  general,  show  cyclogenesis  near  the 
axis  of  a  K  maxima.  As  the  storms  mature,  they  move  away  to  the 
north  and  east,  and  decay  at  a  higher  latitude. 

Figure  11a  which  displays  the  Kg  distribution  for  the  FGGE 
winter,  illustrates  that  care  should  be  taken  in  interpretation  of 
the  results  presented  here.  Note  the  large  eddy  values  off  the 
California  coast.  Recall  that  eddy  values  represent  departures 
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from  the  zonal  mean.  Thus ,  for  example,  if  the  observed  u  component 
was  small  in  an  area  while  u  is  large  for  the  zonal  ring,  the  result¬ 
ing  u*  («u-u)  is  large  negative  which  becomes  very  large  positive 
when  squared.  In  a  sense,  such  a  region  is  an  example  of,  for  lack 
of  a  better  name,  a  K^  "hole."  Large  standard  deviations  of  Kg 
(Fig.  lib)  appear  in  this  region  as  well.  A  review  of  mean  700  mb 
maps  show  ridging  in  December  and  January  with  troughing  in  February 
over  this  area.  (See  Figs.  6  and  7  and  Wagner  (1979).)  This 
variation  in  troughs  and  ridges  produces  the  large  standard  devia¬ 
tions  in  this  region  for  the  winter  season.  This  example  off  the 
California  coast  points  to  the  importance  of  interpreting  K^  plots 
only  after  considering  other  features  in  the  zonal  ring,  inspecting 
K  plots,  and  reviewing  flow  patterns  which  generate  the  energy 
fields. 

Changes  from  winter  to  spring  (MAM)  are  most  evident  in 
spring  K  shown  in  Figure  12.  The  Asian  maximum  has  been  displaced 
east  and  slightly  north  of  its  winter  position  and  has  weakened 
slightly.  It,  and  the  Mediterranean/Arabian  maximum  appear  con¬ 
sistent  with  the  300  mb  K  patterns  of  Lahey  et  al.  (1960).  On  the 
east  coast  of  North  America,  however,  the  winter  maximum  nearly  dis¬ 
appears.  Figure  1  in  Koehler  and  Min  (1983)  shows  that  the  strong 
gradient  of  A  present  in  winter  all  but  disappears  from  the  North 
American  coast  In  spring.  The  barocllnlcity  is  greatly  weakened  and 
cannot  support  strong  winds  aloft.  In  addition,  cyclogenesis  off  the 
U.S.  east  coast  is  lower  for  the  FGGE  spring  than  the  long  term 
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average  as  presented  by  Whittaker  and  Horn  (1981) .  These  three 
things  seem  self  consistent  in  that  weak  baroclinicity  and  weak  winds 
provide  less  favorable  conditions  for  cyclogenesis.  As  with  Figure 
11,  Figure  13  presents  Kg  and  its  standard  deviation. 

Summer  (JJA)  finds  the  continental  east  coast  mairtmn  displaced 
to  their  most  northerly  positions  and  at  their  weakest  Intensity, 
as  shown  in  Figure  14.  Perhaps  the  most  striking  feature  of  summer 
is  the  signature  of  the  Indian  monsoon.  Its  signature  is  present 
in  all  summer  maps.  K  shows  an  isolated  center  of  greater  than 
40  Jkg  1  over  the  Indian  subcontinent  that  is  not  present  in  any 
other  season.  A  standard  deviation  in  K  of  greater  than  30  Jkg  1 
extends  over  the  same  region.  The  greater  than  50  Jkg-1  contour  of 
Kj,  shown  in  Figure  15a  reflects  the  Indian  monsoon.  Both  the 
standard  deviation  of  K  (Fig.  14b)  and  of  Kg  (Fig.  15b)  reflects 
che  time  change  in  the  monsoon.  The  relatively  large  deviations  re¬ 
flect  an  intensifying  of  the  circulation  in  early  summer  followed 
by  a  weakening  in  late  summer.  The  most  likely  features  responsible 
for  the  large  values  are  the  low  level  Somali  jet  and  the  high  level 
Easterly  Jet  which  extends  from  the  Indian  subcontinent  westward 
into  Africa.  As  mentioned  earlier  in  relation  to  the  time-latitude 
variation  of  Kg,  these  features  contribute  to  a  low  latitude  summer 
standing  eddy  maximum  in  the  latitude  bands  of  the  Indian  subcontin¬ 
ent  (See  Fig.  5a). 

Autism  (SOM)  brings  an  Increase  in  K  with  strong  maxima  return¬ 
ing  to  the  North  Pacific  and  North  America  off  the  east  coast  of  the 
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continents.  The  centers  are  still  displaced  north  of  their  mean 
winter  positions,  however.  Figures  16  and  17  show  the  plots  of  K, 

Kg  and  their  standard  deviations. 

This  section  has  shown  some  of  the  differences  in  the  geographi¬ 
cal  distribution  of  kinetic  energy  in  the  FGGE  year.  Noteworthy 
features  include  the  nearly  equally  strong  Asian  and  North  American 
K  maxima  in  winter  and  a  near  total  collapse  of  the  North  American 
maxima  in  spring.  This  weaker  maxima  is  reflected  in  reduced  baro- 
clinicity  in  exact  A  plots  from  Koehler  and  Min  (1983)  and  in  re¬ 
duced  cyclogenesis  in  the  region.  The  most  striking  feature  of 
summer  is  the  strong  signal  of  the  Indian  monsoon. 

d.  Ratios  of  kinetic  to  available  potential  energy 

The  A  calculations  of  Min  and  Horn  (1982)  and  those  for  K  per¬ 
formed  here  allow  ratios  of  various  quantities  to  be  formed,  as  was 
done  by  Peixoto  and  Oort  (1974).  Table  4  contains  monthly  A  from 
Min  and  Horn  (1982)  for  the  approximate  and  exact  A  formulations  and 
the  approximate  A  from  the  five  years  of  data  calculations  performed 
by  Peixoto  and  Oort  (1974).  Peixoto  and  Oort's  values  have  been 
converted  to  specific  units  for  comparison.  Table  2,  presented 
earlier,  contains  the  K  quantities,  and  Table  5  contains  the  ratios 
of  K  and  A  as  calculated  from  the  values  in  Tables  2  and  4. 

Many  of  the  differences  between  the  quantities  of  available 
potential  energy  presented  in  Table  4  have  been  discussed  in  Min 
and  Horn  (1982)  and  Min  (1981),  however,  some  items  should  be  re- 
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Table  4.  Monthly  mean  values  of  available  potential  energy.  Left 
two  columns  are  exact  and  approximate  A  from  Min  (1981) 
and  third  coluan  is  approximate  A  from  Peixoto  and 
Oort  (1974)  (converted  to  specific  units).  Units:  Jkg 


A  h 

«  j _ _ _  " 


Mo. 

Exact 

Approx. 

Peixoto 
&  Oort 

Mo. 

Exact 

& 

Approx. 

Peixoto 
&  Oort 

0 

579.9 

660.1 

686.0 

D 

371.6 

568.2 

584.7 

J 

607.4 

710.0 

721.9 

J 

392.0 

605.8 

607.6 

P 

580.5 

722.1 

717.6 

F 

380.9 

633.1 

614.1 

M 

552.9 

643.5 

674.0 

M 

338.1 

567.5 

584.7 

A 

510.4 

571.5 

536.8 

A 

275.8 

499.3 

461.7 

M 

425.2 

421.7 

400.7 

M 

185.6 

357.2 

329.9 

J 

323.6 

282.2 

269.0 

J 

82.2 

219.6 

204.7 

J 

266.6 

203.8 

202.5 

J 

33.3 

147.1 

141.6 

A 

306.2 

233.8 

218.9 

A 

74.7 

179.0 

163.3 

S 

360.7 

319.7 

320.1 

S 

161.8 

271.0 

259.2 

0 

429.9 

460.8 

457.3 

0 

250.6 

401.5 

381.1 

N 

507.1 

588.7 

597.8 

N 

306.6 

503.9 

504.2 

Mo. 

Exact 

AEtt 

Approx. 

Peixoto 
&  Oort 

Mo. 

Exact 

^Est 

Approx. 

Peixoto 
&  Oort 

D 

64.7 

50.0 

60.3 

0 

139.6 

41.9 

41.1 

J 

73.3 

55.6 

63.7 

J 

142.1 

48.6 

50.5 

F 

64.0 

49.4 

59.8 

F 

135.6 

39.6 

43.3 

M 

54.6 

45.3 

56.2 

M 

160.1 

30.7 

32.7 

A 

71.8 

48.4 

52.3 

A 

162.9 

23.8 

18.1 

M 

70.4 

39.9 

51.7 

M 

169.3 

24.6 

18.9 

J 

56.0 

32.1 

44.1 

J 

185.4 

30.5 

20.5 

J 

49.6 

27.1 

37.5 

J 

183.7 

29.6 

23.7 

A 

56.2 

29.0 

36.9 

A 

175.3 

25.8 

18.7 

S 

64.2 

30.1 

46.3 

S 

134.7 

18.6 

14.5 

0 

67.7 

38.6 

54.0 

0 

111.7 

20.7 

21.6 

N 

70.5 

48.6 

58.6 

N 

130.0 

36.2 

35.5 
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Table  5.  Ratios  of  mean  monthly  kinetic  energy  and  available 

potential  energy  quantities.  The  left  hand  column  pre¬ 
sents  ratios  based  on  exact  A,  while  the  two  right  columns 
show  values  based  on  the  approximate  A  as  calculated  by 
Min  (1981)  for  the  FGGE  year  and  Peixoto  and  Oort  (1974) 
for  a  five  year  period.  For  K-  to  Rz  ratios,  the  left 
column  is  for  the  FGGE  year  ana  the  right  column  is  the 
five  year  average  from  Peixoto  and  Oort. 


K/A 

Peixoto 

VKz 

Peixoto 

Mo. 

Exact 

Approx. 

and  Oort 

Mo.  FGGE 

&  Oort 

D 

.26 

.22 

.24 

D 

1.02 

1.27 

J 

.26 

.22 

.26 

J 

1.08 

1.16 

F 

.27 

.21 

.26 

F 

.98 

1.05 

M 

.26 

.22 

.25 

M 

1.00 

1.09 

A 

.22 

.19 

.25 

A 

1.30 

1.44 

M 

.22 

.22 

.26 

M 

1.75 

2.18 

J 

.21 

.24 

.31 

J 

2.29 

2.78 

J 

.22 

.29 

.38 

J 

2.12 

2.75 

A 

.19 

.25 

.37 

A 

2.20 

2.22 

S 

.19 

.21 

.29 

S 

1.85 

2.11 

0 

.22 

.20 

.24 

0 

1.57 

1.99 

N 

.25 

.21 

.23 

N 

1.63 

1.56 

VAz 

Peixoto 

Mo. 

Exact 

Approx. 

&  Oort 

D 

.55 

.16 

.17 

J 

.55 

.17 

.19 

F 

.52 

.14 

.17 

M 

.64 

.13 

.15 

A 

.85 

.14 

.16 

M 

1.30 

.18 

.21 

J 

2.96 

.29 

.32 

J 

7.11 

.39 

.43 

A 

3.08 

.31 

.34 

S 

1.22 

.18 

.23 

0 

.72 

.15 

.23 
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emphasized  since  they  are  Important  In  explaining  the  differences 
in  ratios  shown  in  Table  5  and  discussed  below.  Very  large  dif¬ 
ferences  occur  in  the  partitioning  of  A  into  zonal  and  eddy  com¬ 
ponents  depending  on  whether  the  exact  or  approximate  formulation 
is  used.  Although  the  annual  cycle  of  total  A  and  Ag  in  the  approx¬ 
imate  form  is  slightly  larger  than  in  the  exact  form,  Ag  in  the 
exact  form  is  nearly  twice  as  large  throughout  the  year,  and  exhibits 
a  slight  summer  maximum  rather  than  a  minimum  as  found  for  the 
approximate  fora.  Ag  is  larger  than  Ag  during  the  summer  months 
while  it  is  always  smaller  in  the  approximate  fora.  The  partitioning 
of  Ag  into  standing  and  transient  components  produces  quite  large  dif¬ 
ferences  between  the  approximate  and  exact  forms.  Agst  is  several 
times  larger  in  the  exact  form  for  all  months.  Note  also  that  while 
Min  (1981)  and  Peixoto  and  Oort  (1974)  agree  quite  well  for  approx¬ 
imate  A,  Ag  and  AEst>  approximate  Agtr  from  Table  4  is  consistently 
lower  by  about  20Z.  This  difference  is  probably  attributable  to 
the  problem  of  interpolating  data  used  in  obtaining  A  to  grid 
points  rather  than  interpolating  the  basic  data.  This  problem  was 

illustrated  in  an  earlier  section  in  relation  to  wind.  Here,  the 

2 

differences  arise  from  T'  quantities  in  Ag  (see  equations  for  A  in 
Min  and  Horn  (1982))  where  T'  is  the  time  average  deviation  from  the 
zonal  average  temperature  on  an  Isobar ic  surface. 

Comparing  the  monthly  ratios  of  K  to  A  in  Table  5  from  Peixoto 
and  Oort  (1974)  with  those  using  approximate  A  determined  here,  one 
notes  the  seasonal  trend  is  similar,  however,  the  ratios  from  this 


study  are  smaller.  This  could  be  due  to  the  different  vertical  do¬ 
mains  of  the  calculations,  different  data  sets  and  possibly  the 
interpolation  techniques.  The  highest  ratios  of  K  to  A  in  both 
studies  occur  in  summer  with  somewhat  lower  ratios  in  winter.  The 
ratio  of  K  to  A  using  A  in  the  exact  form  shows  a  reversal  in  the 
size  of  the  ratios  from  winter  to  summer,  with  larger  ratios  in 
winter  and  smaller  in  summei . 

The  ratios  of  to  found  here  are  in  general  agreement 
with  those  of  Peixoto  and  Oort  (1974).  Except  for  November,  all 
ratios  found  .here  are  lower.  Again,,  this  is  probably  due  primarily 
to  interpolation  differences  in  computing  K£cr  and  the  data  sets 


involved. 
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5.  SPECTRAL  ANALYSIS  RESULTS 

Although  the  focus  of  this  thesis  has  been  on  the  temporal  and 
spscial  distribution  of  kinetic  energy  during  the  FGGE  year,  the 
daily  values  calculated  here  readily  lend  themselves  to  time  series 
analysis.  This  brief  chapter  shows  the  time  series  analyses  of  K, 

A  (both  exact  and  approximate)  and  their  cross  spectra.  The  pro¬ 
cedures  employed  to  perform  this  analysis  were  discussed  in  Chapter 
3;  details  of  data  preparation  prior  to  using  an  EFT  to  determine 
coefficients  for  the  spectral  estimates  are  in  Appendix  B.  Before 
presenting  the  time  series  results,  some  additional  description  of 
the  procedure  employed  is  in  order. 

The  spectral  estimates  presented  below  were  averaged  over  fre¬ 
quency  as  suggested  by  Bendat  and  Plersol  (1971)  to  increase  the 
number  of  degrees  of  freedom  (df) .  Because  a  10Z  cosine  taper 
was  applied  to  each  end  of  the  data  series,  each  spectral  estimate 
has  only  1.74  df  rather  chan  the  normal  2  df.  Three  frequency  bins 
were  chosen  to  form  each  spectral  estimate,  allowing  5.22  df.  A 
Markov  red  noise  spectrum  was  calculated  for  each  series  using  the 
method  discussed  by  Mitchell  et  al.  (1966).  Note  that  this  may  not 
represent  th<*  actual  significance  since  Che  population  to  test  from 
is  not  known.  Nevertheless,  this  provides  a  background  spectrum  for 
comparison.  To  be  consistent  with  the  spectral  estimates,  three  red 

noise  bins  were  averaged  prior  to  multiplying  by  the  appropriate 
2 

(chi-squared)  x  /df  value  for  each  confidence  limit  obtained.  For 
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coherence  significance  levels  of  12  and  52,  Che  three  bin  estimates 
give  values  of  .815  and  .713  respectively.  See  Panofsky  and  Brier 
(1958)  for  this  equation. 

The  resulting  power  spectra  for  K,  Kg  and  Kg  are  shown  in 
Figure  18.  The  solid  line  connects  the  three  frequency  bin  esti¬ 
mates.  The  dashed  lines  are  the  Markov  red  noise  confidence  limits 
for  902,  992  and  99.92.  The  spectrum  is  presented  as  percent  var¬ 
iance  explained  by  each  frequency  band.  All  three  spectra  show  a 
number  of  significant  frequencies,  however,  those  at  23  and  26 
cycles/year  are  especially  striking.  For  K  (Fig.  18a)  the  26 
cycles/year  frequency  exceeds  the  99.92  limit  and  for  Kg  (Fig.  18c), 
both  23  and  26  cycles/year  are  well  above  the  99.92  confidence 
limit.  They  correspond  to  cycles  of  about  16  and  14  days  respec¬ 
tively.  Recall  that  the  visual  inspection  of  Figure  1  suggested 
that  cycles  with  these  periods  were  prominent  in  winter  and  late 
autumn.  They  also  correspond  to  the  periodic  variation  found  by 
Miller  (1974)  during  the  period  January  1964  to  May  1966.  The  re¬ 
sults  shown  in  Figure  18c  indicate  that  those  cycles  appear  even 
when  the  entire  year  is  examined.  The  higher  significant  peak  at 
26  cycles/year  and  also  one  at  32  cycles/year  (near  11.5  days)  occur 
in  Kg  while  these  frequencies  in  the 
Despite  this,  cross  spectral  analysis  between  K^  and  Kg  shows  a 
coherence  significant  at  or  above  12  for  these  frequencies  as  shown 
in  Figure  19.  (The  dashed  lines  are  the  12  and  52  significance 
levels.)  Phase  angle  difference  (A8),  not  shown,  indicates  that  at 


Kg  series  are  not  significant. 
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Figure  18. 


Power  variance  spec t run  for  K.  Expressed  as  percent 
variance  explained.  Dashed  lines  are  confidence 


limits  of  90,  99  and  99. 9X.  a)  K.  b)  K_.  c)  Kg. 
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figure  19.  Coherence  from  cross  spectral  analysis  of  Kg  and  Kg. 
Dashed  lines  are  5 Z  and  1Z  significance  levels. 

the  14  day  cycle  Kg  leads  Kg  by  135*  (5*25  days)  and  at  11.5  days 
Kg  leads  Kg  by  152°  (4.86  days).  This  agrees  with  our  understand-* 
ing  of  the  energy  cycle  which  shows  the  flow  of  energy  from  Kg  to 

V 

McGuirlc  and  Reiter  (1976)  noted  strong  periodicities  at  about 
24  days  in  several  energy  quantities  from  a  time  series  covering 
nearly  ten  consecutive  years.  They  noted  the  periodicity  was 
strongest  in  the  winter  months  and  tended  to  disappear  in  summer. 
With  regard  to  their  findings,  the  FGGE  year  would  appear  to  be 
anomolous.  No  strong  periodicity  around  24  days  was  found.  The 
coarseness  of  resolution  at  lower  frequencies  and  possible  time 
series  errors  introduced  by  Interpolating  for  missing  values  may 
have  contributed  to  not  observing  their  cycle.  Another  explanation 
could  be  that  the  FGGE  year  was  quite  different  from  other  years 
with  regard  to  this  particular  periodicity. 
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Because  dally  A  values  were  available  for  the  period  of  this 
study,  these  time  series  were  also  analyzed.  Spectral  variance 
estimates  were  prepared  for  exact  and  approximate  A  and  are  pre¬ 
sented  in  Appendix  C.  Coherence  estimates  between  the  two  A  forms 
and  K  quantities  are  also  included.  The  plots  are  presented  with 
interpretation  left  to  the  reader. 


6.  CONCLUSIONS 


Because  kinetic  energy  is  so  directly  related  to  atmospheric 
motion,  its  seasonal  and  geographical  distribution  is  of  particular 
interest  to  the  meteorologist.  This  study  has  focused  on  twice 
daily  kinetic  energy  calculated  for  the  Northern  Hemisphere  using 
the  level  Ilia  FGGE  year  (December,  1978-November ,  1979)  data  set. 
Computations  were  performed  in  isentroplc  (6)  coordinates.  Total 
kinetic  energy  (K)  was  divided  into  its  zonal  (Kg)  and  eddy  (Kg) 
parts,  and  the  results  presented  as  a  time  series  of  twice  daily 
values.  Kg  was  further  divided  into  its  standing  and  transient  corn* 
ponents.  The  geographical  distribution  of  vertically  Integrated 
grid  point  values  of  K  and  Kg  were  calculated  along  with  their 
standard  deviations.  By  using  exact  and  approximate  available  poten¬ 
tial  energy  (A)  values  from  Min  and  Horn  (1982) ,  ratios  of  K  to  A 
were  evaluated.  Finally,  time  series  of  K  and  A  were  subjected  to 
spectral  and  cross  spectral  analysis.  The  following  major  results 
were  obtained: 

(1)  As  expected,  time  series  of  twice  daily  values  show  quite 
strong  annual  cycles.  Kg  shows  the  weakest  annual  cycle.  Kg  ex¬ 
hibits  a  more  pronounced  annual  cycle.  Visual  inspection  reveals 
that  during  winter  and  again  in  late  autumn  a  strong  pulsing  occurs 
at  a  period  of  12  to  16  days. 

(2)  Comparison  of  transient  eddy  K  with  vslues  from  Pelxoto  and 
Oort  (1974)  showed  tbs  FGGE  values  to  be  lower  by  about  27%.  While 


year  to  year  differences  exist,  a  smaller  vertical  domain  and 
systematically  lower  estimates  resulting  from  a  different  inter¬ 
polation  method  probably  are  the  main  contributors  to  the  lower 
values . 

(3)  The  plot  of  standing  eddy  K  as  a  function  of  time  and 
latitude  exhibited  three  relative  maxima  during  the  year.  The  first 
was  a  high  latitude  maximum  associated  with  persistent  east-west 
trough  and  ridge  features.  An  extreme  case  was  observed  in  the  5* 
latitude  bands  between  70-80°  N  during  February.  It  resulted  from  a 
very  strong,  persistent  wave  number  two  pattern.  The  second  maxi¬ 
mum  was  present  throughout  the  FGGE  year  and  was  closely  aligned 
with  the  mean  monthly  jet  stream  position.  The  third  maximum 
occurred  in  summer  at  about  15®N.  This  feature  appeared  to  be  a 
monsoon  signal. 

(4)  The  time  and  latitude  plot  of  transient  eddy  K  appeared 
about  5  to  10°  north  of  the  second  standing  Kg  maximum,  and  in 
general,  the  latitude  of  greatest  cyclogenetic  activity.  Noteable 
exceptions  occurred  during  February  and  May  when  the  maximum  fre¬ 
quency  of  cyclogenesis  was  displaced  north  of  the  transient  Kg  maxi¬ 
mum.  A  comparison  of  transient  Kg  with  transient  eddy  available 
potential  energy  from  Min  and  Born  (1982)  showed  the  transient  Kg 
maximum  north  of  the  exact  transient  Ag  maxima  while  it  was  south 
for  approximate  transient  Ag. 

(3)  Geographical  plots  of  grid  point  K  showed  the  Asian  and 
North  American  east  coast  maxima  to  be  nearly  equal  strength  during 


the  FGGE  winter  (DJF) ,  while  the  long  term  mean  Indicates  the  Asian 
maximum  is  the  stronger.  As  expected,  the  upper  level  maxima  are 
closely  associated  with  regions  of  cyclogenesis  in  winter  as  well 
as  the  other  seasons.  A  very  dramatic  change  occurred  during  spring 
(MAM):  the  North  American  maximum  nearly  disappears.  A  weak  grad¬ 
ient  of  A  shown  in  Koehler  and  Min  (1983)  and  decreased  cyclogenesis 
compared  with  the  20  year  average  from  Whittaker  and  Horn  (1981) 
appear  consistent;  with  the  greatly  reduced  North  American  K  maximum. 
The  most  striking  feature  of  simmer  is  the  signature  of  the  Indian 
monsoon. 

(6)  Ratios  of  K  to  A  using  exact  A  show  the  ratio  highest  in 
winter  and  lowest  in  simmer.  This  is  the  reverse  of  that  found 
using  approximate  A.  A  lower  ratio  of  K  to  approximate  A  between 
Pelxoto  and  Oort  (1974)  and  this  study  seems  related  to  domain  and 
interpolation  differences  as  well  as  yearly  variability. 

(7)  Spectral  analysis  produced  strong  periodicity  in  Kg  at 
14  and  16  days.  This  confirmed  earlier  results  obtained  from 
visual  inspection.  Cross  spectral  analysis  of  Kg  and  Kg  revealed 
highly  significant  coherence  at  14  days  with  Kg  leading  Kg  by  5.25 
days.  The  strong  24  day  cycle  found  by  McGulrk  and  Reiter  (1976) 
in  nearly  ten  years  of  data  was  not  evident  in  any  of  the  K  or  A 
time  series.  This  would  seem  to  indicate  the  FGGE  year  was  quite 
different. 

This  study  has  raised  some  interesting  points  for  future  in¬ 
vestigation.  The  relationship  between  K,  A  and  cyclones  should  be 
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explored  further.  A  study  involving  a  limited  region  of  strong 
FGGE  cyclone  activity  might  reveal  more  concrete  links  between 
cyclones  and  their  energy  fields.  Comparisons  between  the  winter 
and  spring  season  off  the  east  coast  of  North  America  could  add  in¬ 
sight  into  the  large  seasonal  change  in  this  region.  Additionally, 
time  series  analysis  of  twice  daily  energy  quantities  at  selected 
grid  points,  in  latitude  bands  or  summed  over  limited  regions  might 
reveal  more  about  the  periodicities  discussed  here  on  a  hemispheric 


scale 
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APPENDIX  A 


This  appendix  briefly  discusses  the  procedure  used  to  calculate 
the  area  of  grid  squares,  and  those  used  to  perform  the  vertical 
integration. 

The  area  of  grid  squares  is  determined  by  integration  of  the 
equation 

ds  ■  a^  cos  <p  d$d\  (Al) 


where  a  is  the  mean  radius  of  the  earth,  <J>  is  latitude  and  A  is 
longitude.  Integration  of  Equation  Al  gives 
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4^  and  <j>s  are  the  north  and  south  latitude  limits  of  a  grid  square. 
For  calculations  here,  ^  and  are  5*  apart.  A$  and  AA  are  also 
5*.  The  term  sin  $/A$  constitutes  the  map  factor.  This  method  of 
area  calculation  provides  an  exact  determination  of  the  grid  square 
area  because  it  results  from  the  direct  integration  of  Equation  Al. 
This  method  is  more  accurate  than  one  method  which  uses  the  average 
of  the  cosine  of  the  north  and  south  latitudes  (cos  $),  and  another 
which  uses  the  cosine  of  the  average  of  the  north  and  south  lati¬ 
tude  (cos  $). 

Vertical  integration  of  the  kinetic  energy  equations  for  each 
grid  square  is  performed  using  the  trapesoldal  rule.  The  value  of 
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u,  v,  and  p  on  two  isentropic  surfaces  separated  by  A0  -  10K  is 
shown  in  Figure  AI,  The  upper  and  lower  values  of  the  component 

Q 

winds  are  averaged  (u  -  (uu  +  11^/2,  etc.)  while  Ap  is  given  by 
Ap  ■  PL  "  Pu*  This  gives  layer  average  wind  speeds  and  a  mass 
increment.  This  procedure  is  performed  over  the  domain  from  as 
low  as  220K  to  370K.  Summation  over  the  actual  number  of  vertical 
layers  yields  units  of  energy  per  unit  area  (Jm~2)  for  the  grid 
square.  This  value  can  then  be  divided  by  the  mass  of  the  atmos¬ 
phere  in  the  domain  of  integration  to  yield  specific  units  (Jkg-1). 


0+  A0 


Uu.  VU»  Pu 


*  U®,  V® 


6 - UL,VL,PL 

Figure  Al.  Schematic  of  interpolation  of  layer  mean  quantities 
used  in  vertical  integration.  Subscripts  refer 
to  quantities  at  6  surface  data  levels.  Ap  is  pressure 
difference  for  layer,  and  u°  and  v®  are  layer  average 
winds. 
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APPENDIX  B 

This  appendix  briefly  discusses  Che  procedures  employed  prior 
to  spectral  analysis  of  kinecic  energy  and  available  potential 
energy  time  series. 

A  fast  Fourier  transform  (FFT)  routine  available  at  the 
Madison  Academic  Computing  Center  (see  MACC,  1973)  and  based  on  the 
Cooley  and  Tukey  (1965)  algorithm  was  employed.  Prior  to  using  the 
FFT,  certain  requirements  must  be  satisfied.  For  example,  there 
must  be  no  missing  values  in  the  series.  Consequently,  the  time 
series  gaps  (730  periods  -  699  available  ”  31  gaps)  were  filled  by 
linear  interpolation  of  values  from  the  nearest  available  adjacent 
time  periods  on  each  side  of  the  missing  period.  This  resulted  in 
a  time  series  of  length  730.  Next,  factoring  requirements  for  the 
FFT  routine  must  be  met.  Since  728  was  a  convenient  number  of  data 
periods  that  met  the  factoring  requirements,  the  last  two  time 
periods  were  arbitrarily  dropped.  These  remaining  728  periods 
were  subjected  to  a  binomial  filter  of  weights  .25/. 50/. 25  (Hanning 
filter)  to  reduce  diurnal  variability  which  was  found  to  be  large 
by  Koehler  and  Min  (1983)  for  values  of  Ac.  Although  diurnal 
changes  were  small  for  K  quantities,  it  was  decided  to  process 
all  the  time  series  in  the  same  manner. 

Because  the  FGGE  period  is  of  finite  length,  certain  adjust¬ 
ments  must  still  be  made.  For  computational  purposes,  an  FFT  as¬ 
sumes  a  time  series  can  be  wrapped  back  on  itself  in  a  differen- 
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Ciably  veil  behaved  manner.  To  enaure  Che  beginning  and  end  of  a 
series  meet  correctly,  the  series  mean  is  removed  and  the  first 
and  last  l/10th  of  the  series  is  subjected  to  a  cosine  taper  con¬ 
sisting  of  a  left-  and  right-half  cosine  bell  (see  Bingham  et  al. , 
1967).  This  process  reduces  leakage:  a  problem  whereby  a  pure 
wave  of  given  frequency  is  represented  by  a  number  of  frequencies 
after  transformation.  Finally,  the  series  is  standardized.  This 
ensures  that  the  variance  sums  to  one  and  allows  spectral  esti¬ 
mates  at  each  frequency  to  be  readily  interpreted  in  terms  of  Z 
variance  explained.  For  additional  information  on  the  above  proce¬ 
dures,  see  Bergland  (1969),  and  references  therein,  Bendat  and 
Piersol  (1971)  and  Bloomfield  (1976).  Two  specific  applications 
to  meteorology  are  found  in  Madden  and  Julian  (1971)  and  Julian 
(1971) . 
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appendix  c 

This  appendix  shows  the  spectral  variance  plots  for  exact  and 
approximate  available  potential  energy.  In  addition,  coherence 
plota  are  presented  to  show  the  results  of  cross  spectral  analysis. 
In  each  pair  of  coherence  plots,  exact  A  quantities  appear  in  the 
top  diagram  while  approximate  A  quantities  are  shown  below  for 
comparison.  As  discussed  in  Chapter  5,  the  dashed  lines  at  .815  and 
.713  are  the  significant  levela  of  1Z  and  5Z.  Interpretation  of 
these  figures  is  left  to  the  reader. 
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Figure  C3.  Coherence  from  cross  spectral  analysis  of  A  and  K 
Dashed  lines  are  1Z  and  5Z  significant  levels, 
a)  Exact  A  and  K.  b)  Approximate  A  and  K. 
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